Continuous spin reorientation in antiferromagnetic films 
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We study anisotropic antiferromagnetic one-layer films with dipolar and nearest-neighbor ex- 
change interactions. We obtain a unified phase diagram as a function of effective uniaxial D e and 
quadrupolar C anisotropy constants. We study in some detail how spins reorient continuously below 
a temperature T s as T and D e vary. 
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Considerable attention has been devoted to the mag- 
netic properties of ultrathin magnetic films in the last 
years. An interesting feature of magnetic films is dis- 
continous spin reorientation (DSR), i.e., thermally driven 
switching between perpendicular and in-plane spin align- 
ment at a temperature T r below the ordering tempera- 
ture To. DSR has been studied both experimentally and 
theoretically. It has been established that DSR depends 
on the competition between dipolar interactions and uni- 
axial anisotropy often found in films. Continuous spin 
reorientation (CSR) is also very interesting. A thermally 
driven CSR transition was first observed experimentally 
in bulk systems at some temperature T s well below To[l|. 
Below T s , all spins rotate continuously as a whole as T is 
varied in these systems. After these experiments Horner 
and VarmaQ proposed an early phenomenological model 
in which higher anisotropes, which compete with the uni- 
axial anisotropy, were required for obtaining CSR. More 
recently, CSR has been observed in ferromagnetic thin 
films |3j- Some nonhomogeneus multilayer models have 
been proposed to explain CSR0. 

The aim of this paper is to report numerical results 
on CSR for one-layer antiferromagnetic films. It is im- 
portant to note that CSR is always associated with a SR 
phase defined by its own broken symmetries 0, 0. For 
instance, the order parameter m may be perpendicular 
(9 = 0) to the film plane in the T s < T < T range, 
and tilt away (0 < 9 < n) from the easy magnetization 
axis below T s , thus breaking additional symmetries 0. 
To our knowledge, such SR magnetic phase has not been 
observed in numerical simulations of antiferromagnetic 
films. 

We consider a system of classical unit spins {S ; } in a 
square lattice with hamiltonian H = Hj + Hd + H a where 
Hj, Hd, and H a are for short range exchange, long range 
dipolar, and anisotropy interactions, respectively. We use 
periodic boundary conditions. The exchange and dipo- 
lar energy between two antiparallel out of plane nearest 
neighbors spins is J < and — Sd, respectively. Fur- 
thermore, there is site uniaxial —D(Sf) 2 plus fourfold 



— C[(Sf) + (Sf) 4 ] anisotropy energies. Similar models 
have been studied for C = 0, and some DSR between de- 
generate in-plane and out-of-plane states have been found 
as T, J and D vary0. For ferromagnetic films, on the 
other hand, Jensen et al. have considered C ^ but 
obtained phase diagrams only for T = 0. 

Ground state configurations obtained from Monte 
Carlo simulations are exhibited in Fig. 1. States desig- 
nated with a c, can be defined @ by Sf = r? cos 9, Sf = 
rf sin sin <j>, Sf — rf sin 9 cos (f> where rf — (—l) y ^\ 
T ? = and rf = (-1)^)+^). On the other 

hand, states designated with a s are defined by rf = 
rf = rf = (— l^W+^M. A suitable order parameter for 
both c and s is m a = N~ x £\ Sf rf . 

We calculate energies for these configurations as in 
Ref. [6] , and find a surface anisotropy energy A that be- 
haves as an easy axis anisotropy. We find that A = 2sd 
for s states and A = — 1.23ed + 2J for c states. This sug- 
gests we can define an effective anisotropy as D e = D + A 
and obtain a unified phase diagram for both s and c states 
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FIG. 1: Magnetic phase diagram for dipolar systems in their 
ground states for J < 0. s and c states correspond to exchange 
dominated systems with J < — 1.6l£d and dipolar dominated 
systems with — 1.61e<j < J < 0, respectively. Full and dashed 
thick lines stand for first- and second-order transitions, re- 
spectively. SR stands for the spin reorientation phase, in 
which < 9 < it/2. 
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FIG. 2: Tilt angle 9 versus T/Tq for a pure dipolar system 
(J = 0) of 32 x 32 x 1 spins with C = —Ed- The system is 
ordered for T < To. Curves correspond from top to bottom 
to D e /C = 0.9,0.8,0.7,0.6,0.5, and 0.4 respectively. Similar 
plots have been obtained for J < 0. In the inset, 9 versus 
~D e /C for T — is shown. Solid line corresponds to tan9 = 
\f D e / (C — D e ). x (o) stand for systems of 32 x 32 spins with 
J = and variable D (D = and variable J), respectively. . 

for T = 0, as shown in Fig. 1. s configurations give a 
lower energy for J < — 1.6l£d while c states are more 



favorable for — 1.61e<2 < J < 0. Interesting experimen- 
tal realizations of the latter condition could be found in 
Ref. [9]. 

In both cases we find a z-collinear phase (9 = 0) if both 
D e > and D e > C are fulfilled. In-plane configurations 
(9 = 7r/2) are prefered for D e < C. More interestingly, 
we obtain a spin reorientation phase for C < D e < 
in which 9 covers the < 9 < ir/2 range. Minimiza- 
tion of the total energy gives tan9 = -J D e / (C — D e ) 
and therefore spins rotate continuously from 9 = to 
7r/2 as D e varies from to C, as shown in the inset 
of Fig. 2. Symbols in the same inset correspond to nu- 
merical data obtained by cooling from the paramagnetic 
phase to T « Ed, J for different values of D e /C. 

We have explored temperature driven CSR in the SR 
phase by MC simulations for both exchange and dipo- 
lar dominated systems. For that purpose we calculate 
the order parameter (m x ,m y ,m z ) and the energy as a 
function of T. We find two different regions (see figs. 2 
and 3). Upon cooling below To, in-plane configurations 
appear for D e /C > 0.65, and spins rotate towards the 
z axis below a second order transition at T s < To. On 
the other hand, for D e /C < 0.65 spins point out of the 
plane below T , and rotate towards the xy plane as T de- 
creases below T s . Finally, we find, that the ratio T s /Tq, 
as in d = 3 systems @, seems to depend mainly on D e /C 
and not on J or Sd- 
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FIG. 3: Phases of pure dipolar (J = 0) films for C = —Ed 
obtained from MC simulations, o (x) stand for systems of 
16 x 16 (32 x 32) spins respectively. Transition tempera- 
tures have been obtained from peaks observed in the specific 
heat. Systems have been cooled in AT = — 0.01ed/fcs steps 
of 4 x 10 5 MC sweeps each. Similar phase diagrams have been 
obtained for J < 0. 



